Quantum read-out and fast initialization of nuclear spin qubits with electric currents 
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Nuclear spin qubits have the longest coherence times in the solid state, but their quantum read- 
out and initialization is a great challenge. We present a theory for the interaction of an electric 
current with the nuclear spins of donor impurities in semiconductors. The theory yields a sensitivity 
criterion for quantum detection of nuclear spin states using electrically detected magnetic resonance, 
as well as an all electrical method for fast nuclear spin qubit initialization. 
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Running a current over a localized electron spin pro- 
vides an extremely sensitive method to detect electron 
spin resonance [l[ . This technique, commonly referred to 
as "electrically detected magnetic resonance" (EDMR) is 
now experiencing a revival in the context of solid state 
quantum computing Several interesting recent re- 
sults [1 -12 point to the possibility of using EDMR for 
quantum read-out of individual donor nuclear spin states 
in silicon. The hope is that EDMR will be the "silicon 
alternative" to optically detected magnetic resonance, a 
method that can detect single nuclear spin states of de- 
fects in diamond [i~3| . Nuclear spins in silicon are excep- 
tional quantum memory units, with coherence times Ti n 
over a few seconds hj, |l_5| . While single spin detection of 
donor electrons in silicon was demonstrated recently fl6| , 
the corresponding read-out of individual nuclear spins re- 
mains a great open problem. 

Nuclear spin read-out using EDMR Q is based on the 
interaction between donor electron (S) and nuclear (J) 
spin in an external magnetic field B: 



H, 



AS- 1, 



(1) 



where ui e = g e ji e B/H and u> n = g n [i n B/h are the Zee- 
man frequencies for electron and nuclear spins, respec- 
tively. The hyperfine interaction A couples the donor 
electron to its nuclear spin. From the approximation 
1-L en ~ (uj e + AI Z )S Z , we see that the electron spin res- 
onance frequency will depend on the nuclear spin state. 
The EDMR experiment consists in applying microwaves 
at a fixed frequency wesr, and measuring an electri- 
cal current as a function of the applied magnetic field 
B. Thus, the current will show a single peak either 
at = S,(o;esr — A/2) / (g e [i e ) when the donor nu- 
clear spin is aligned along the B field, or at B + = 
^(^esr + ^/2)/(<?eMe) when the donor nucleus is aligned 
against the field. It is therefore believed that the detec- 
tion of electron spin resonance of a single donor impurity 
using ultra sensitive EDMR methods will allow the read- 
out of its nuclear spin state, or even of 29 Si nuclear spins 
nearby the donor. 

However, it is not clear that this is possible. The 
problem is that the time it takes to detect EDMR of 
a single donor spin tEDMR must be shorter than nu- 



clear spin- flip time T\ n , and as we show below the elec- 
tric current will induce nuclear spin-flips. Moreover, 
read-out time must be longer than nuclear spin coher- 
ence time for the nuclear spin wave function to 
collapse into one of the outcome states. Thus, the 
EDMR read-out scheme will be limited by a contrast 
of exp [— Max(T2„,iEDMR)/7i n ]. The best possible sce- 
nario occurs when tEDMR < T2 n , leading to an optimal 
read-out contrast of exp (— T^n/Tin). 

In this letter, we describe our theory for the interaction 
of nuclear spin qubits with a "classical current" , i.e. a 
current formed by an incoherent ensemble of conduction 
electrons, the one used in common electronic devices such 
as EDMR. 

Consider a donor impurity interacting with an elec- 
tron gas [3-12|. Scattering events preserve the total spin 
of conduction and localized electrons, leading to the in- 
teraction Hamiltonian 17 
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Here c\ a ,Ck a are creation/annihilation operators for a 
conduction electron with momentum k and spin a =t 
or 4-, S± = S x i iS y are spin raising/lowering operators 
for the donor electron, and J is their "bare" exchange 
interaction. 

Exchange scattering will tend to equilibrate the local- 
ized spin with the electron gas. If the localized spin is not 
in thermal equilibrium, its expectation value (S a ) along 
any direction a = x, y, z will decay exponentially in time 
according to the rate, 
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where v is the energy density of states at the Fermi level. 
The effective exchange interaction J e g is obtained after 
summing and renormalizing all orders of perturbation 
theory in Ji>, resulting in an universal expression that 
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FIG. 1. (color online) This figure demonstrates how a donor 
impurity implanted in a silicon accumulation field effect tran- 
sistor (inset) can be tuned into and out of the Kondo regime 
using a top gate voltage. The plot shows | In (T/Tx) for 
T = 0.1 K and T = 1 K for a single donor impurity fixed at a 
distance R from the transistor interface (model assumptions 
and parameters are the same as Fig. 7 of @|). The length 
scale zo models the two-dimensional electron gas width, that 
is controllable electrically over the range zo ~ 30 - 100 A. 
Hence a donor can be tuned from R/zo ~ 3 (strong coupling 
regime) to R/zo w 5 (weak coupling regime) by increasing the 
top gate voltage. Recently, a similar electrical tuning of the 
donor Kondo temperature was demonstrated in a side-gated 
silicon device 12011. 



depends only on the Kondo temperature Tk [171 1 1 81 ] : 
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Calculations of the Kondo temperature for realistic pa- 
rameters of a silicon transistor are presented in Figure [T] 
Equation © is valid for T > T K . When T < T K , 
the conduction electron spins will screen out the donor 
electron spin, forming a Kondo singlet. This will add 
non-exponential decay to (S a ), with a rate r e ~ ksTx /h 

At T > Tk, (S a ) decays exponentially in time, towards 
its thermal equilibrium value. This leads to the follow- 
ing effective model for donor electron plus nuclear spin 
evolution subject to the conduction electron environment 



dp 

dt 
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Here p is the density matrix describing donor electron 
plus nuclear spin, r e is given by Eq. ([3|), and p e = 
— tanh [huj e / (2kBT)] is the equilibrium donor electron 
spin polarization. With the usual definitions (S) = 



Tr{pS}, (I) = Tr{pl}, and (SI) = Tr (pSI) (a ma- 
trix formed by the outer product between S and I) we 
can put Eq. ([5]) in a more convenient form that we call 
the generalized Bloch equation for the donor: 

(S) = uo e z x(S)-A(SxI)- T e ((S) - ^z) (6a) 

( j) = -u n z x (/) + A (S x 7) , (6b) 
(SI) = u) e z x (SI) + Lu n (SI) x z 
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with £ = e Q( g 7 e Q ,e / 3e 7 the Lcvi-Civita tensor. 

We can solve Equations (|6a | - (|6"c"]) analytically by mov- 
ing to a reference frame where the electron spin rotates 
with frequency lo £ and the nuclear spin rotates with — co n ; 
within order (A/ui e ) 2 we get 
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where (I') is the nuclear spin operator in the rotating 
frame, and B = (uj e + ui n ). The decay rate of Eq. ([Taj) 
is the new nuclear spin dccoherence rate induced by the 
current. 
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In the limit T e <C A, we get \/T-2n ~T e /2 showing that 
the nuclear spin "follows" the electron. However, in the 
opposite regime of T e 3> A we get a decoherence rate 
that is inversely proportional to T e : l/T 2n « A 2 /T e , 
signaling a transition to a motional narrowing regime. 
Here the electron spin fluctuates so fast that its effect on 
the nuclear spin's phase is averaged out. 

Similarly, the decay rate of Eq. (|7b[) is the new nuclear 
spin-flip rate, 
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As T e increases, 1/Ti n peaks at T e — B w ui e . The origin 
of this effect is a transfer of population from the electron 
spin to the nucleus assisted by the "flip-flop" coherence 
(S-I+). It is necessary for the electron spin fluctuation 
frequency T e to match the flip-flop coherence B for the 
transfer to occur. The net result is a transfer of electron 
spin polarization p e into the nuclear spin. 

This result suggests a method for qubit initialization 
via interaction with the electron gas. The idea is to 
electrically tune the Kondo temperature so that T e gets 
close to ui e . This can be achieved when In (T/Tk) ~ 
^371-/(2^1) or T « 10 T K for \p e \ = 1, and we call 
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this the "strong coupling regime" (r e ps cj e ). Fig- 
ure [T] describes how to tune into this regime with a 
top gated device. In the strong coupling regime we get 
Max(l/T ln ) ps A 2 /(4w e ). For a phosphorus donor in sil- 
icon, A = 120 MHz and cu e = 28 GHz for B = 1 T, 
leading to qubit initialization times T\ n = 1 /is. These 
electrical initialization times are several orders of mag- 
nitude faster than other methods previously considered 
for donor nuclear spins in silicon. In 22( , white light was 
used to polarize donor nuclei within 150 s, while |23j used 
a laser tuned to an exciton transition to initialize within 
0.5 s. 

We remark that the nuclear spin-flip process described 
by Eq. (|7b[) results from an "indirect" interaction of the 
donor nucleus with conduction electrons, mediated by 
the donor electron spin through hypcrfine coupling A. 
There also exists an additional "direct" interaction be- 
tween the donor nuclear spin and the conduction elec- 
trons: It gives rise to the Korringa relaxation mechanism 
24j, that is usually quite weak for semiconductors. We 
calculated the Korringa rate for the phosphorus donor 
in the silicon transistor shown in Fig. [TJ We obtained 
(1 /T\ n ) Korringa < 10~ 3 s" 1 for R/z > 2 at low tem- 
perature (T < 10 K). Hence, at times shorter than a few 
minutes we may neglect this "direct" interaction between 
the donor nucleus and the electron gas, and the donor nu- 
clear spin polarization will remain close to the electron's, 
Pn ~ Pe = — tanh (huj e /2kBT). However, we note that 
the Korringa rate will eventually drive the nuclear spin 
back to its thermal equilibrium state, with small polar- 
ization p n = + tanh {hjj n /2k,BT). 

We now describe nuclear spin read-out with EDMR. 
The donor spin is irradiated with a resonant microwave 
field of frequency wesr and amplitude B± , inducing elec- 
tron spin precession with frequency f2j_ = g e ^ e B±/h. 
Only electron spins satisfying the equal population or 
"saturation" condition fl± > T e can be detected [l[, 
and optimal sensitivity is achieved at the threshold for 
saturation, fl± = T e Q. In terms of the Kondo tem- 
perature, optimal EDMR takes place for In (T/Tk) = 
y[37r/ (2\p e \)](B/B±), a regime that requires low Kondo 
temperature. We call this limit the "weak coupling" 
regime. For B = 1 T and B± = 0.3 G, weak coupling 
requires \n(T/T K ) ~ 500. 

We generalized Eqs. (|6a[) - ([6c)) for the case of resonant 
excitation; when Q±_ > T e and Q±_ <C A, the nuclear 
spin relaxation rates are given by Eqs. (|Taj) - (|7b[) with 
r e — > Qj_ . Hence we have T ln ~ 2{uj e /A) 2 /f2j_ an d 
T2n ~ 2/17 j_. From this result we can deduce the maxi- 
mum possible contrast for read-out of nuclear spin states: 
Nuclear spin-flips during read-out will lead to contrast 
equal to exp (— TmjTin) = exp [— (A/ui e ) 2 ], that is quite 
close to 1 when uj e 3> A. We recall that this optimum 
contrast occurs only when £edmr, the time it takes to 
detect EDMR of a single donor, is less than T2 n . 

A more likely scenario is that £edmr will be quite long. 



Assuming the longest acceptable £edmr is equal to Ti n , 
we now derive a sensitivity condition for EDMR detection 
of nuclear spins. Such a condition will depend critically 
on the value of the EDMR current amplitude (AJ)edmr, 
that is given by the current difference when only one 
donor is on and off resonance with a microwave field. A 
key observation is that (AJ)edmr must be larger than 
the shot noise (A/) sno t accumulated during T\ n : 
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where I is the average device current and N(T\ n ) = 
ITi n /e is the number of electrons that passed through 
the device during T\ n . Using the value of T\ n computed 
above we get the sensitivity criterion for nuclear spin 
read-out: 
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We note that an implicit assumption of Eq. (|lll) is that 
A > 1/T 2 *ei i- e -j the separation of the two EDMR peaks 
is larger than their lincwidth. For 7 = 1 /iA, B± = 0.3 G 
and B = 1 T we get (A///) E dmr > 3 x 10~ 6 

A recent theory of the spin-dependent scattering mech- 
anism of EDMR predicted (A//7) ED mr = 6x 10~ 6 x 
tanh(fiw e /2/c B T) x (B/VT) x (Bj_/0.3G), for a single 
donor implanted in a transistor of area 0.1 (/mi) 2 . At 
B ~ 1 T and T = 1 K this is right at the sensitivity 
limit; however, much better sensitivity can be obtained 
by going to higher B fields. 

There exists other methods for EDMR detection that 
can potentially yield much higher (A///)edmrj most no- 
tably the ones based on optical excitation and recombi- 
nation P, UH, 0] ■ Electron- hole pairs are excited by light 
and the photocurrent is monitored as a function of mag- 
netic field. In this case additional relaxation channels 
such as capture of an extra donor electron with rate r c , 
with corresponding ionization of the extra electron with 
rate Ti can be present. It is straightforward to generalize 
our theory to these cases; we obtain the same expressions 
as above with the substitution T e — > (T e +T t +T C ) and 
A 2 -> (A 2 ) = A^iKT^Tc) (effective hyperfine coupling 
gets reduced by capture, because a singlet produces zero 
hypcrfine field on the nuclear spin). With these substi- 
tutions, Eqs. (|8j)- (fTTj) can be applied to recombination 
based EDMR. 

In a typical transistor, the gate voltage can be used 
to control the overlap between the electron gas and the 
donor impurity wavefunction. Thus, the Kondo tempera- 
ture can be controlled electrically from Tk = (when the 
electron gas density is zero, i.e., the transistor is turned 
off) all the way to large Tk- This is demonstrated for a 
top gated device [l2|, [25| in Figure Q] Hence, we argue 
that it is possible to switch donors electrically into and 
out of three distinct modes: In the "quantum evolution" 



4 



mode, donors have no overlap with the electron gas, im- 
plying Tjc = and T± n , T2 n are longer than a few seconds 
[la] . For larger gate voltage, overlap is small and the 
donor is in the weak coupling regime, allowing quantum 
read-out. For even larger gate voltages, the donor will be 
in the strong coupling regime enabling fast initialization. 

In conclusion, we presented a theory for the interac- 
tion of an electric current with donor nuclear spin qubits. 
We showed that electric tuning of donors away from the 
Kondo regime allows coherent evolution and high con- 
trast quantum read-out using electrically detected mag- 
netic resonance; in contrast, electric tuning close to the 
Kondo regime enables fast initialization of nuclear spin 
qubits. We derived a general sensitivity criterion for the 
read-out of nuclear spin states using EDMR. 
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